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A transformation of supersymmetric quantum mechanics for coupled channels is presented, which 
allows the introduction of up to TV degenerate bound states without altering the remaining spec- 
trum of the Hamiltonian. Phase equivalence of the Hamiltonian can be restored by two successive 
supersymmetric transformations at the same energy. The method is successfully applied to the "^Si- 
^Di coupled channels of the nucleon-nucleon system and a set of Moscow-type potentials is thus 
. . , generated. 

Q . PACS numbers: 03.65.Nk, 13.75.Cs 

O ■ 

o ■ 
(N : 

=3 . 



< 



> 



O 

o 



X 



I. INTRODUCTION 



The formalism of supersymmetric quantum mechanics (SQM), introduced by Witten [y in 1981, provides an elegant 
way to construct a hierarchy of Hamiltonians with well defined relations between their spectra ||l|-^. Specifically, 
00 it allows the elimination or introduction of bound states |^ without altering the remaining part of the spectrum. 
, Transformations of SQM, which we will refer to as supersymmetric (SUSY) transformations, are specific Darboux 
transformations of the Schrodinger equation and are related to the factorization method . Consequently, there 
is a close connection to inverse scattering theory PJlO[|. 

The SUSY-transformations of the radial Schrodinger equation have been studied in detail ||^,^ . In particular emphasis 
, was given to mathematical aspects as well as applications to specific physics phenomena. An important application is 
' the relationship of equivalent effective interactions between composite particle systems. In such systems, because of 



00 the necessary suppression of the internal degrees of freedom, the resulting interactions are ambiguous [gT 12 1. There 
is always a set of effective interactions leading to the same scattering phase shifts which sustain, however, a different 
number of additional unphysical bound states as a consequence of simulating differently the Pauli principle. 
Baye pointed out that SQM is an elegant way to construct such phase equivalent potentials. Since then, more 
■ general transformations of the radial Schrodinger equation for uncoupled channels have been worked out |^J^ which, 
I , apart from a simple removal and an addition of a ground state, allow arbitrary modifications of the bound state 
' Q spectrum jl^ ]. 

^ The concept of SUSY-transformations has been extended to coupled-channel systems by Amado et al. but 
their transformations do not allow the construction of phase-equivalent potentials after the removal of a bound state 
Recently, Sparenberg and Baye ||l^ addressed this problem and presented SUSY-transformations for coupled- 
channel systems which in addition to the removal of bound states allow the construction of phase-equivalent potentials. 
However, the status of SUSY-transformations for coupled-channel systems is still incomplete because the process of 
^ introducing a bound state has not been considered so far. In this article we present a SUSY-transformation for the 
- - introduction of N degenerate bound states in a system of N coupled channels. 



II. SUPERSYMMETRIC TRANSFORMATION 



We consider a system of A'^ coupled channels which is described by the Schrodinger equation 



where e = k"^ — 2mE/h^ and Uo{r) = 2mVQ{r)/h'^ . The potential Vb(r) is an A^ x matrix that may include the 
centrifugal barrier and thresholds which may be different in each channel. For simplicity, we assume that the mass 
m 



is equal for all channels (for unequal masses, see Ref. |15|). The wave function ^'o(e,r) = {ipQ{e,r), ■ ■ ■ ,ipQ{e,r)) 
is an N X N matrix, where each column vector ipKe, r) is a solution of Eq. (|^). 
The SUSY transformations are based on the factorization of the Hamiltonian 

Ho = A+A-+e, (2) 
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where the energy e is smaller than or equal to the ground state eq of Hq and 



^ ±— + W(r) . 
ar 

The superpotential W{r) is an N x N matrix satisfying the differential equation 



dW 
dr 



= U, 







The supersymmetric partner Hamiltonian Hi is given by 



A-A+ +e with Ui{r) = Uo{r) 



dr 



(3) 



(4) 



(5) 



The solutions of the corresponding Schrodinger equation, at any energy e, are directly given in terms of 'I'o(e,7') by 
4*1 (e,r) = A^^'o(e,r). The above equations are formally equivalent to those for uncoupled channels for which the 
SUSY transformations are known in closed form ^-|^. In the present work we consider a SUSY transformation for 
coupled-channel systems generated via the ansatz 



Wir) = r)%'{e, r) + i^^lil, r)) ' A 

divl/-i(e-i)(vl/t(e-i)] A 



(6) 



where we have used the hermiticity of the first term of Eq. (||) . The constant N x N matrix A must be symmetric 

and should be chosen in such a way that W{r) has no singularities for r > 0. The first term of Eq. (^) is a 

straightforward extension of the simplest expression for the superpotential in uncoupled systems. 

The superpotential (|^) can also be put in the simplest form where 4'o is now the most general solution of the 

Schrodinger equation JTzt . The associated transformation depends on the choice of 4'o(e, t"), the integration bound a, 

and the constant matrix A. In the present work we are concerned with choices that lead to the addition of a bound 

state. 

For an energy e < eo, it turns out that the regular solution of Hq together with the bound a = oo is the appropriate 
choice. Then the transformed wave function 'I'i(e,r) at the energy e = e, is simply given by 



«'i(e,r) = 1'o(e,r) 



1 - 



dt^-, 



A 



(7) 



The regular solutions 5*0, at energies below the lowest ground state energy of Hq, exhibit a diverging behavior 
for r —>■ oo. Hence, all columns of \E'i(e, r) are exponentially vanishing solutions of Hi at asymptotic distances. 
Furthermore ^'i(e, r) is bounded (for < r < oo) if 



det 1 



dt* 



(8) 



which is the case when the matrix A is chosen negative semidefinite. This guarantees also that Ui{r) is bounded 
except for r = 0. 

The behavior of ^'i(e, r) at r requires more attention as it depends on the coupled-channel system considered. 
Specifically, if the singularity of the potential at the origin is of the form v{v + l)/r^ with the same v for all channels, 
one obtains the behavior lim^^o r)r~^ = const ^ 0. Hence, for v >1 the function 4'i(e, r) vanishes at the origin 
and with an appropriate choice of A, satisfying (8|), each column vector of \E'i(e, r) corresponds to a bound state of Hi. 
Thus the transformation associated with Eq. ) enables us to add N degenerate bound states to the spectrum of Hq 
at the energy e. More precisely, analytical examples show that the degeneracy aX e — e seems to depend on the rank 
of the matrix A For v — the function \E'i(e, r) does not correspond to bound states. The situation becomes 
even more intriguing, when the potential exhibits different singularities in the coupled channels, i.e. z/^, i = 1, . . . , TV. 
It can be shown that if the coupling vanishes near the origin, the ith component of each column vector of 5'i(e, r) 
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behaves as r"^* there and leads, with an appropriate choice of A, to the addition of up to N degenerate bound states 
to the spectrum of Hq. 

For a nonvanishing couphng near the origin and different the behavior of 5'i(e, r) near the origin cannot be given 
in closed form. Therefore, definite conclusions about the number of bound states that can be added are difficult to 
draw. However, the numerical examples discussed below demonstrate that the transformation (|^) can also lead, in 
the general case, to the addition of up to N degenerate bound states. 



III. PHASE-EQUIVALENCE 

The SUSY transformation generated by W{r) of Eq. (||) modifies, similarly to the single channel case, the S-matrix 
To compensate this, another transformation is required which is associated with a superpotential that has the 
opposite sign at asymptotic distances fl^ . Therefore we first perform a transformation generated by the simple 
superpotential (first term of Eq. (^)) using asymptotically vanishing solutions 770 of the Schrodinger equation 
with Ho [lim,_^oo ?7o(e, f ) exp(\/l r) = A, where ^4 is an x matrix]. The existence of the transformed potential 
requires a non- vanishing determinant of ?7o(e, r) for all finite r- values. This can always be satisfied because e is assumed 
to be smaller than the deepest bound state of Hq. The first SUSY transformation modifies the S-matrix but does 
not change the number of bound states. In a second step we transform Hi via the transformation mediated by the 
superpotential of Eq. (^) . 

A direct verification shows that {rio^ {£,r)y is a regular solution of Hi. We use this solution to construct the 
superpotential of the second transformation, following the principle of the preceding section. The potential U2 
resulting from the two successive transformations, as well as the corresponding solution ^'2, can be expressed in terms 
of the solutions of the initial equation only. They read 

U2{r) = t/o(r) - 2^ [x2{e,rWo{e,r)} (9) 

and 

*2(e, r) = -(e - e)^'o(e, r) - X2(e, r)W[j]Q{e, r), *o(e, r)] , (10) 



where 



X2{e,r) = ?7o(e,r)A 



1- / dtr,l{e,t)rj„{e,t)A 



(11) 



corresponds to the added bound state(s) of H2 aX e — e. The quantity A is a constant N x N matrix. 
This procedure transforms a regular solution \l/o into a regular solution \E'2 as long as the singularity at the origin 
is sufficiently strong {i^ > 2). In this case the new potential U2 is phase equivalent to Uq [Eq. (^o|) shows that both 
^0 and ^2 have the same asymptotic behaviour, up to a normalization factor, which confirms phase equivalence], 
and sustains (besides the bound states of Hq) up to N additional bound states at e. Similarly to the single SUSY 
transformation, discussed in the previous section, the degeneracy of the introduced bound states at e can be controlled 
by the number of non-vanishing eigenvalues of the constant A 
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IV. EXAMPLES 



As a demonstration we consider an £ ^ 2 two-channel system without threshold and a potential of Gaufi form, 
T/y = exp(-rVi?2) + 5ij6h^/2mr^, with R ^ 2fm, 2mc^ = 938.9185MeV and the depths f/"' = -lOOMeV, 

^2^2^ = — 60MeV, and y/2'' — ^2? ~ — 56.56 MeV. This potential does not sustain a bound state. A pair of degenerate 
bound states at ii^ = —200 MeV has been introduced using the twice iterated SUSY transformation. The new potentials 
are shown in Fig. 1 together with the corresponding wave functions of the two added bound states. It is seen that 
there is a clear spatial separation between the two degenerate states which is also reflected in the potential matrix. 
Choosing A of rank one results in the introduction of only one bound state at e. This leads to the potential and 
bound state wave functions displayed in Fig. 2. It must be emphasized that the potentials of Fig. 1 and Fig. 2 are 
phase equivalent. 

The introduction of bound states via the twice iterated SUSY transformation (p|-pT|) is also possible for systems where 
channels with different quantum numbers of orbital angular momentum are coupled. Of particular interest is the 
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coupling of an s-channel (£i = 0) and a d-channel {I2 = 2) when the associated potential has only in the d-channel 
an 1/r^-singularity at the origin. In such a system the strength of the singularity suffices for the introduction of 
one additional bound state via the SUSY transformation The transformed potential has no r^'^ singularity 

at the origin which, taking into account the centrifugal term of the d-channel explicitly, implies an attractive 
contribution to the residual potential. We have verified this special case numerically assuming the previous two-channel 
system with slight modifications, i.e. i?i = 0,i?2 = 2 and Vij ~ V^p exp(— r^/i?^) + (52j6?i^/2mr^. This system has a 
bound state at i? = —29.45 MeV. The introduction of an additional bound state at E — 200 MeV leads to the potential 
and bound state wave function displayed in Fig. 3. The transformed potential confirms the previous discussion and 
its phase equivalence has been verified numerically. On the contrary, for a A matrix of rank 2, we have verified that 
the final potential is not phase equivalent to Uq. 

As a realistic example we consider the '^Si-'^Di partial waves of the nucleon-nucleon (NN) system which are coupled 
by a tensor potential that sustains a single physical bound state at Ed — —2.22 MeV (the deuteron). Within the quark 
model the NN force is an effective interaction of a composite particle system and may giv e rise to Pauli-forbidden 
states which are simulated by deep unphysical bound states in local potential models [ pl| . Neudatchin et al. [T^ ] 
were the first to use the concept of forbidden states for the NN potential. A more detailed study led to the well 
known Moscow potential which reproduces the scattering data up to 400 MeV Elimination of the additional 

unphysical bound state via the twice iterated SUSY transformation of Sparenberg and Baye leads to a phase 
equivalent potential having a close similarity to the standard NN interactions derived in the meson-exchange picture. 
Qualitatively the same result has been obtained by Leeb et al. [po| by applying a numerical SUSY transformation to 
an early version of the Moscow potential. 

In the present work, we apply the procedure described by Eqs. (|9|-pTl) in the opposite direction. Starting from the 
Reid soft core potential (RSC) ||2^ for the ■^Si-^Di channel, we would like to generate a series of NN potentials of 
Moscow-type where an additional bound state has been introduced at different energies e. In Fig. 4 we show a series of 
phase equivalent NN-potentials obtained by this procedure. Variation of the matrix A results in modifications of the 
r-dependence while the energy of the added bound state is reflected in the depth of the potentials. The repulsive core 
of the effective potential V22{r) is a consequence of the repulsive core of the RSC potential. It must be emphasized 
that the obtained potentials are not exactly of the same nature as the Moscow potential since the central part of their 
D-channel component has an r^^ attractive core. 



V. CONCLUSIONS 



In summary, in this work we have presented a closed form expression for a SUSY transformation which enables us to 
introduce single or degenerate bound states in a coupled channel system. The novel transformation works for a general 
hermitean A^-channel system with and without thresholds. The method has been successfully tested for the case of 
two coupled partial waves in the NN system where a whole set of phase equivalent NN potentials of Moscow-type has 
been generated. The method has also been successfully applied to two-channel systems with potentials of Gauf5 form. 
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FIGURE CAPTIONS 



Figure 1 

Introduction of two degenerate bound states at E = — 200MeV into a two-channel system with Gaufiian potential 
(see text) at ^ = 2 using An = —5 • 10*, A22 = —5 • 10^, A12 = A21 = and ?7o(e, r) associated with An — A12 — 1 and 
A21 = —^22 — 0.5. The upper graph shows the matrix elements of the potential; Vii(r) (solid hne), V22('') (dashed 
line) and Vi2{r) = V2i(r) (dotted line). The two lower graphs show the components of the wave functions (solid line 
for the first channel, dashed line for the second channel) of the two added bound states. 

Figure 2 

The same system as shown in Fig. 1 but entering only one bound state a.t E — — 200MeV using An — —5 • 10'', 
A12 ~ A21 — A22 = and the same tjq. For the notation see Fig. 1 

Figure 3 

Introduction of a bound state a.t E = — 200MeV into a two-channel system of coupled s- and d-partial waves with 
Gaufiian potential (see text) using An = —5 • 10*, A12 = A21 = A22 = and ryo(e, r) associated with An = A12 — 1, 
A21 = —A22 = 0.5. The upper graph shows the matrix elements of the potential; Vii(r) (solid hne), V22{r) (dashed 
line) and Vi2{r) — V2i{r) (dotted line). The graph at bottom shows the components of the wave function (solid line 
for the first channel, dashed line for the second channel) of the added bound state. 

Figure 4 

NN potentials for the ^Si-'^Di-channel, phase equivalent to the RSC generated by the twice iterated SUSY transfor- 
mation (||). The unphysical bound states are at E — — lOOMeV (long dashed line), E = — 200MeV (dashed line), 
E = -400 MeV (dotted line) and E = -800 MeV (sohd line). The constant matrix A was fixed to A22 = -5 • 10*^, 
All = A12 = A21 = and rjo{e,r) associated with An = A12 = 1 and ^21 = —A22 = 0.5. 
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